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Конечно-элементное моделирование пьезоэлектрического устройства 
накопления энергии на основе кантеливера* 


А. Н. Соловьёв, Ле Ван Зыонг 


Рассматривается моделирование пьезоэлектрического генератора, который является элементом устройства 
накопления энергии. Пьезоэлект, рический генератор представляет собой консольно-закреплённую пластину, 
на которую наклеены пьезоэлектрические элементы и присоединена инерционная масса. Исследуются два 
варианта возбуждения колебаний пластины. В первом — закреплённая сторона совершает вертикальные 
гармонические колебания на заданной частоте и с заданной амплитудой, во втором к ЭТОЙ СТ ороне приложе- 
на сила, изменяющаяся по гармоническому закону. В качестве математической модели УСТ! ройст, ва рассмат- 
ривается трехмерная краевая задача линейной теории электроупругости для составного упругого и ПЬзо- 
электрического тела. Решение краевой задачи проводится методом конечных элементов в пакете АМ5у5. 
При численном решении в качестве пьезоэлектрического материала взята пьезокерамика ПКР-7М, а в каче- 
стве материала пластины рассмотрены стеклопластик, дюраль, сталь. В качестве материала инерционной 
массы используется алюминий. Численно исследованы два случая колебания на резонансной частоте и ко- 
лебания в низкочастотной области на частоте значительно меньшей частоты первого резонанса. Проведено 
исследование зависимости резонансной частоты устройст! ва от толщины пластины для различных материа- 
лов и от величины массы инерционного элемента. Результаты представлены в виде Г| рафиков, ПОЗВОЛЯЮЩИХ 
найти резонансную частоту для определённых размеров. Исследованы зависимости выходного потенциала 
на свободных электродах пьезоэлементов на резонансных частотах и в низкочаст отной области от этих же 
параметров, эти результаты представлены также в виде Г! рафиков, что позволяет конструктору выбрать рацио- 
нальные размеры элементов и сочетание материалов для достижения наибольшей эффект, ивности уст! ройст! ва. 
Ключевые слова: МКЭ, накопление энергии, пьезоэлектрик, оптимизация, кантилевера. 


Введение. Современная малогабаритная бытовая техника, телефоны сотовой связи, беспровод- 
ные сенсорные системы для наблюдения и диагностики технического состояния различных объек- 
тов и т. п., потребляют весьма незначительную энергию. Для питания такого рода устройств ши- 
роко применяются автономные источники электрической энергии, которые в ряде случаев могут 
быть заменены пьезоэлектрическими преобразователями энергии совершающими механических 
колебания. 

По сравнению с другими методами накопления энергии, такими как электромагнитные [4] 
и электростатические [5], пьезоэлектрические накопители энергии на базе преобразования коле- 
баний окружающей среды привлекают большое внимание в последнее время из-за их простой 
структуры, прямого преобразования энергии колебаний в электрическую энергию с высоким 
уровнем напряжения [6, 7]. 

Наиболее известны пьезоэлектрические преобразователи энергии двух конфигураций: 
осевого и кантилеверного типа. Они используются как для создания пьезогенераторов электриче- 
ской энергии, входящих в состав автономных источников питания, так и для высокоточных ли- 
нейных двигателей или пьезоэлектрических актюаторов. Пьезоэлектрическое устройство накоп- 
ления энергии кантилеверного типа работает эффективнее, когда частота возбуждения находится 
в непосредственной близости от основной резонансной частоты электромеханической системы. 

Некоторые упрощенные аналитические модели для пьезоэлектрического накопителя энер- 
гии кантилеверного типа представлены в литературе. Например, в [8] разработана аналитическая 
модель с распределёнными параметрами. Вычисление энергии пьезоэлектрических структур при- 





* Работа выполнена при финансовой поддержке РФФИ (гранты № 13-01-00196-а, 13-01-00943-а). 
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ведено в [9]. 5бодапо и др. использовали метод Релея — Ритца и разложения по собственным 
формам колебаний для описания устройства накопления энергии [10]. ЕйигК и Тптап в [11] по- 
строили для одномерной модели устройства накопления энергии решение в аналитической фор- 
ме. Цао и 5о4дапо [10, 12] модифицировали модель 5одапо, чтобы иметь возможность эффективно 
изучать системы накопления энергии и определять их оптимальные параметры. 

В настоящей работе изучено влияние размеров и условия закрепления устройства накоп- 
ления энергии при возбуждении колебаний на резонансной частоте или вблизи резонансной ча- 
стоты, а также в низкочастотной области вдали от резонансных частот. Решение проблемы осно- 
вано на конечно-элементном моделировании в пакете АМ$У$. 

На рис. 1 представлено устройство накопления энергии с использованием биморфной 
конфигурации. Структура кантилеверного типа пьезоэлектрического устройства накопления энер- 
гии, состоит из металлической пластины (подложки) 2, консольно закреплённой на основании 4, с 
двумя приклеенными пьезоэлементами 1 в зоне максимальных механических напряжений, возни- 
кающих при изгибе. С целью повышения коэффициента преобразования за счёт снижения резо- 
нансной частоты на верхней поверхности подложки может быть дополнительно размещен инер- 
ционный элемент 3. Устройства накопления энергии с консольно-закреплёнными изгибными эле- 
ментами спроектированы для преобразования внешней энергии окружающей системы в низкоча- 
стотном диапазоне. 
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Рис. 1. Структура пьезоэлектрического устройства накопления энергии на основе кантеливера: а — общая структура, 
6 — первая модель, в — вторая модель. 1 — пьезоэлемент; 2 — подложка; 3 — элемент инерционный; 4 — основание 


При анализе собственных колебаний пьезоэлектрического устройства накопления энергии 
предполагается выполнение следующих электрических и механических граничных условий. Плос- 
кие поверхности пьезопластин покрыты тонкими электродами, представляющими собой эквипо- 
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тенциальные поверхности. Вся внешняя поверхность пьезопреобразователя свободна от механи- 
ческих напряжений. 

1. Континуальные постановки задач электроупругости. Рассмотрим некоторый пьезопре- 
образователь ©, представленный набором областей ©, =@„; К=1,2,...,№М,; Л=К со свой- 
ствами пьезоэлектрических материалов и набором областей © =; Т=Ь2,....М; 
1=М, +т со свойствами упругих материалов. Будем считать, что физико-механические процес- 


сы, происходящие в средах О и можно адекватно описать в рамках теорий пьезоэлектри- 


ет ! 
чества (электроупругости) и упругости. 
Для пьезоэлектрических сред © ‚= О предположим, что выполняются следующие поле- 


вые уравнения и определяющие соотношения [13]: 


р„и+а,р,а-У.в=Ё,;У.09=0, (1) 
сс". (Е+ВуЕ)-е/ Е; О-+с.Ю=е, - (#+6,2)+э'.Е, (2) 
= = (Уи+ Уи) /2; Е=-УФ, (3) 


где р(х) — плотность материала; и(х, #) — вектор-функция перемещений; в — тензор механи- 
ческих напряжений; # — вектор плотности массовых сил; В — трёхмерный вектор индукции 
электрического поля; с’ — тензор четвертого ранга упругих модулей, измеренных при по- 


стоянном электрическом поле; е, — тензор пьезомодулей третьего ранга; = — тензор де- 
формаций; Е — трёхмерный вектор напряжённости электрического поля; ф(х,Ё) — функция 
электрического потенциала; эх — тензор второго ранга диэлектрических проницаемостей, 


измеренных при постоянной деформации; а,,Ву, с. — неотрицательные коэффициенты 


демпфирования, а остальные обозначения стандартны для теории электроупругости, за ис- 
ключением дополнительного индекса 7 указывающего на принадлежность к среде О, с номе- 


ром 1. 
Для сред ©, =®„„ с чисто упругими свойствами будем учитывать только механические 


поля, для которых примем аналогичные (1)—(3) полевые уравнения и определяющие соотношения 
в пренебрежении электрическими полями и эффектами пьезоэлектрической связности. 

Система дифференциальных уравнений (1)—(3) дополняется механическими и электриче- 
скими граничными условиями и среди них отметим условие на электроде $, который связан с 
электрическим устройством накопления энергии. 


[2,%=1, (4) 


где Г — ток в цепи, который в случае свободного электрода равен нулю. 
2. Конечноэлементное моделирование. В конечноэлементном пакете АМЗ\У$ осуществлено 
построение модели пьезоэлектрического устройства накопления энергии на основе кантеливера с 
подложкой, выполненной из стеклопластика (дюраль или сталь) размером /[хВхН. Геометриче- 
ские размеры приведены на рис. 1. Пьезоэлементы (ПЭ) состоят из двух одинаковых пьезопла- 
стин, поляризованных по толщине, имеют соответствующие размеры /,х ВьхН,. 

В настоящей работе рассматривается пьезоэлектрическое устройство накопления энергии 
для двух моделей с разными условиями закрепления и механического воздействия. 
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Первая модель (рис. 1, 6) совершает малые колебания в подвижной системе координат, 
связанной с поверхностью, к которой прикрепляется левая сторона пластины. Вертикальные ко- 
лебания подвижной системы определяются соотношением: 


у = у, - ехр(-/. (217) -Ё), (5) 
где у — амплитудное колебание (у = 0,1 мм), Ё— частота вынужденных колебаний. 


Во второй модели основание устройства может свободно перемещаться по вертикали и 
закреплено в горизонтальном направлении (рис. 1, в). Колебания вызываются действием силы РЁ: 


Е=Е, ехр(-/. (2п/)-#), (6) 


где А — амплитудная внешняя сила (А, = 1 Н), /— частота вынужденных колебаний. 

На размеры конструкции (рис. 1) приняты следующие ограничения: размеры пластины 
[хВ = 110х10 мм; толщина пластины #Н изменяется в пределах 0,05-2 мм; пьезоэлементы имеют 
размеры [»хВьхН, = 54х6х0,5 мм; величины [м и Ви постоянны и равны /„ = 8 мм; В» = 22 мм; 
значения /л1 = [5 = 1 мм. 

Ниже приводятся константы материалов, используемые в расчетах: пьезоэлемент 1 — 
пьезокерамика ПКР-7М; подложка 2 — стеклопластик, дюраль, сталь; инерционный элемент 3 — 
алюминий. 











Таблица 1 

Материальные константы пьезокерамики ПКР-7М 
СЕ,ГПа | СЁ,ГПа | СЁ,ГПа | СЁ,ГПа | СЁ,ГПа | е„,Кл/м? | ез,Кл/м? | е,,Кл/м? | ЕЕ, = /& 
109 67 64 97 23 —13,1 27,1 20,1 1990 1740 
































= =8, 85х10? Ф/м, плотность р = 8,0х103 кг/м?, добротность О = 10. 

Упругие свойства изотропных материалов задаются с помощью модуля Юнга Ё и коэффи- 
циента Пуассона У. 

Рассматривается преобразователь с подложкой из различных материалов: дюраль, сталь и 
стеклопластик. Значение плотности дюраля выше значения плотности стеклопластики, но ниже 
значения плотности стали. 

















Таблица 2 
Характеристики конструкционных материалов 
Стеклопластик Дюраль Сталь Алюминий 
Плотность, р (кг/м?) 2500 2700 7800 2700 
Модуль Юнга, Е (ГПа) 70 74 210 68 
Коэффициент Пуассона, \ 0,25 0,34 0,3 0,36 

















В данной работе для двух моделей, исследуются значения собственной резонансной ча- 
стоты #.и амплитудные значения электрического потенциала Уна верхнем свободном электроде, 
в зависимости от материала и значения толщина подложки Р при отсутствии и наличии инерци- 
онной массы, а также в зависимости от ее значения М. 

3. Результаты численных расчётов. 

3.1. Значения толщины подложки А изменяются в пределах 0,05 - 2 мм, без инерци- 
онной массы. В этом случае (инерционная масса отсутствует) на рис. 2 представлены результа- 
ты зависимости значений собственной частоты #. и амплитудные значения электрического потен- 
циала Уна верхнем свободном электроде при колебаниях на резонансных частотах от материала 
и значения толщины подложки А, изменяющейся в интервале 0,05 - 2 мм. 
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Рис. 2. Зависимости собственной частоты от материала и значения толщина подложки р без инерционной массы: 

а, в — для первой модели; 6, г— для второй модели. 1 — стеклопластик; 2 — дюраль; 3 — сталь 


На рис. 2, а, 6 представлены зависимости собственной частоты от материала и значения 
толщина подложки /. 

Из рис. 2, а, б видно, что для обеих моделей, значения собственной частоты возрастает с 
ростом толщины. Значения собственной частоты для первой модели меньше (рис. 2, а), чем для 
второй модели (рис. 2, 6). При этом собственные частоты при материале подложки — сталь выше 
частот для подложки из стеклопластика, но ниже частот для подложки из дюраля. 

На рис. 2, в, гпредставлены зависимости амплитудных значений электрического потенци- 
ала на верхнем свободном электроде от материала и значения толщина подложки А, при колеба- 
нии на первой собственной частоте. 

Из рис. 2, в видно, что для первой модели значение электрического потенциала на верх- 
нем свободном электроде возрастает с ростом толщины В. И наоборот, из рис. 2, г видно, для 
второй модели значение электрического потенциала на верхнем свободном электроде убывает с 
ростом толщины подложки 1. 

3.2. Значения толщины подложки АЙ изменяются в пределах 0,05 -- 2 мм и инерцион- 
ная масса равна 5 г. В этом случае на рис. 3 представлены результаты зависимости значения 
собственной частоты и амплитудные значения электрического потенциала на верхнем свободном 
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электроде при колебаниях на резонансных частотах от материала и значения толщины подложки 
В, изменяющейся в интервале 0,05 - 2 мм. 
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Рис. 3. Зависимости собственной частоты от материала и значения толщина подложки А, инерционная масса равна 5 г: 
а, в— для первой модели; 6, г— для второй модели. 1 — стеклопластик; 2 — дюраль; 3 — сталь 
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На рис. 3, а, 6 представлены зависимости собственной частоты от материала и значения 
толщины подложки А. 

Из рис. 3, а, б видно, что для обеих моделей, значения собственной частоты возрастает с 
ростом толщины. Значение собственной частоты для первой модели (рис. 3, а) меньше, чем для 
второй модели (рис. 3, 6). При этом собственные частоты при материале подложки — дюраль 
выше частот для подложки из стеклопластика, но ниже частот для подложки из стали. 

На рис. 3, в, гпредставлены зависимости амплитудных значений электрического потенци- 
ала на верхнем свободном электроде от материала и значения толщина подложки А, при частоте 
колебаний равной значению первой собственной частоте устройств. 

Из рис. 3, в видно, что в первой модели, значение электрического потенциала на верхнем 
свободном электроде возрастает с ростом толщины 1. Наоборот, из рис. 3, гвидно, что во второй 
модели, значение электрического потенциала на верхнем свободном электроде убывает с ростом 
толщины 1. 

3.3. Значения инерционной массы М изменяются в пределах 0 - 15 г и толщина под- 
ложки Л равна 1,2 мм. В этом случае результаты зависимости значения собственной частоты и 
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амплитудные значения электрического потенциала на верхнем свободном электроде при колеба- 
ниях на резонансных частотах от материала и значения массы инерционного элемента М, изме- 
няющейся в интервале 0 - 15 г, представлены на рис. 4. 
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Рис. 4. Зависимости собственной частоты от материала и значения массы инерционного элемента М, А = 1,2 мм: 
а, в— для первой модели; 6, г— для второй модели; 1 — стеклопластик; 2 — дюраль; 3 — сталь 


На рис. 4, а, 6 представлены зависимости собственной частоты от материала и значения 
массы инерционного элемента М. 

Из рис. 4, а, 6 видно, что для двух случаев, значение собственной частоты убывает с ро- 
стом толщины. Значение собственной частоты для первой модели меньше (рис. 4, а), чем для 
второй (рис. 4, 6). 

На рис. 4, в, г для первой и второй моделей соответственно представлены зависимости 
амплитудных значений электрического потенциала на верхнем свободном электроде от материа- 
ла и значения массы инерционного элемента М при колебаниях на первой собственной частоте 
устройств. 

Из рис. 4, в видно, что в первой модели значение электрического потенциала на верхнем 
свободном электроде убывает с ростом массы инерционного элемента М. И наоборот, из рис. 4, г 
видно, что во второй модели значение электрического потенциала на верхнем свободном элек- 
троде возрастает с ростом массы инерционного элемента М. 
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3.4. Значения электрического потенциала на верхнем свободном электроде при коле- 
баниях на нерезонансной частоте 1.0 Гц. В этом случае рассматривается задача возбуждения 
колебаний на нерезонансной частоте 10 Гц. При этом значения электрического потенциала на 
верхнем свободном электроде, в зависимости от материала и значения толщины подложки 1, при 
отсутствии и наличии инерционной массы, а также в зависимости от её значения М представлены 
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Рис. 5. Зависимости амплитудных значений электрического потенциала на верхнем свободном электроде от материала и 
значения толщины подложки 1: а, в, д— для первой модели; 6, г, е— для второй модели: 1 — стеклопластик; 
2 — дюраль; 3 — сталь 
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На рис. 5, а, 6, ви гпредставлены зависимости амплитудных значений электрического по- 
тенциала на верхнем свободном электроде от материала и значения толщины подложки [, при 
колебаниях на нерезонансной частоте 10 Гц. Рис. 5, а, 6 описывают случай, когда инерционная 
масса отсутствует, а на рис. 5, в, гпредставлены результаты для инерционной массы равной 5 г. 
На рис. 5, д, е представлены зависимости амплитудных значений электрического потенциала на 
верхнем свободном электроде при нерезонансной частоте колебаний 10 Гц от значения массы 
инерционного элемента М, изменяющейся в интервале 0 - 15 гв случаях, когда толщина подлож- 
ки р равна 1,2 мм. 

Из рис. 5 видно, что для двух моделей значение электрического потенциала на верхнем 
свободном электроде, при котором значение частоты колебаний равно 10 Гц убывает с ростом 
толщины А (рис. 5, а, 6, ви Г) и возрастает с ростом значения массы инерционного элемента М 
(рис. 5, д, 6). 

Заключение. В работе рассматривается трехмерная конечноэлементная модель в пакете АМ5У$ 
устройства накопления энергии на основе кантеливера для двух способов возбуждения колеба- 
ний: в первом осуществляется кинематическое воздействие, во втором — силовое. Активным эле- 
ментом устройства накопления энергии является пьезоэлектрическая керамика. Исследована эф- 
фективность устройства в зависимости от его геометрических характеристик и материалов в 
условиях некоторых ограничений на размеры устройства. Проведённые в работе расчёты позво- 
ляют выбрать рациональные размеры элементов и материалы для достижения наибольшей эф- 
фективности устройства. При этом проанализированы два случая. В первом — устройство накоп- 
ления энергии работает на первой собственной частоте изгибных колебаний, а во втором — ча- 
стота колебаний фиксирована и достаточно мала. В обоих случаях в работе сформулированы ре- 
комендации конструктору для выбора материала подложки, её толщины и массы инерционного 
элемента. Построенные зависимости позволяют сконструировать устройство с заданной частотой 
собственных колебаний, которая в реальных условиях может определяться внешним воздействием. 
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ЕУМТТЕ-ЕЕЕМЕМТ МОРЕН-ТМС ОЕ РТЕРОЕТЕСТВТС ЕМЕВСУ $ТОВАСЕ ОЕ\МТСЕ 
ВА$ЕО ОМ САМТТЕЕМЕВ^ 


А. М. Зоомуем, Ее Уап Оцопа 


тишайпд а р/егое/есётс депегайюог и/Р/СН [5 ап еетепЕ о! те епегду 5югаде аемгсе 15 сопх/аегеа, А риегое/есёис 
депегаюг арреаг5 а сап еуег-тоиге4 рае оп и/Р/сВ {те риегое!есёс еетепё аге д/ие, апа пе тега! та55 
5 аНасреа. Тио орйоп$ оЁ те ре иБгаНоп ехсйаНоп аге пиуезНдаеа. [пе ПгзЕ сазе, Те Нхеа $!4е ехесшех 
Веауе Вагптогс тойопз аЁ а дМеп тедиепсу апа ий" {те аезгей атрШМиае; т пе 5есопа — а ваптопк Гогсе 15 
арр/иеа ю 65 $е. А гее-атепзюпа! Боипаагу-иа/ие ргоМет о! {Те /пеаг теогу оЁ еесгое/а$Всйу Гог {пе 
сотрозйе ес апа регое/есёис Боду 15 сопеге4 аз а тайетайса! тодЕ! о! {пе деуксе. Тпе Боипаагу ргоБ- 
[ет [5 сое Егоидй {те Япйе-еетепЕ раскаде АМ5У$. Ипаег {ве питейса/ зоиНоп, РКЕ-7М риегосегат/с 15 
{акеп аз а регоеес тс тейепа| ирегеаз ПЬего/ас$, дига!ит/п, 5! аге сопИАегеа аз а рае таепа/. Аиптйпит 
[5 изе4 аз а тейепа! о! пе тег#а/ таз$. Ти саез аге зи ей сотриайопеи\. [п {те ИтзЕ сазе, пе Нхеа $!ае 
ехесие$ уегйса! вагтогис ифгайоп$ аЁ а омеп ‘едиепсу апа и! а ргезЕЁ атрШиае, т {пе 5есопа сазе, те 
Югсе иагутд т а ваптоп ГазРоп 15 арр!еа ю 5 ясе. Тре аеисе гезопапсе йедиепсу дерепаепсе ироп {те 
р/йе исКпес$ Гог уайоиз таепа/5 апа ироп ве уа/ие оЁ {Те тега! еетепЕ таз$ /5 $НиЕц. ТПе гези!5 аге рге- 
5етеа дгар!иса/у {ПаЁ а!ои5 Рпатд {ве гезопапЕ тедиепсу Юг сепат хез. Тре оиёриЕ роепйа/ дерепаепсе$ аЁ 
Ше Гее еесодез оГ те р/егое!есёс еетепЁ оп {те гезопапЕ тедиепс!е5 апа т те /ои/-Редиепсу гедоп оп 
{е зате рагатеег$ аге пуезНдайе4. Трезе геи аге а/зо ргезегтеа одгарса/у ёпаЁ а/ои5 а Че$дпег {0 5е/есЕ 
геазопа Ме Геа{иге 7е5 апа пе тих оЁ тейепа/5 ю орйтЁе ве авуке. 

Кеуигогабв: ГЕМ, епегду зюгаде, р/егоеесс, орйт!гавоп, сапеуег 
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